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We study the dynamic pliase diagram of a spin model associated with the Number Partitioning 
Problem, as a function of temperature and of the fraction K/N of spins allowed to flip simultaneously. 
The case K — 1 reproduces the activated behavior of Bouchaud's trap model, whereas the opposite 
limit K = N can be mapped onto the entropic trap model proposed by Barrat and Mezard. In 
the intermediate case 1 <^ K <^ N , the dynamics corresponds to a modified version of the Barrat 
and Mezard model, which includes a slow (rather than instantaneous) decorrelation at each step. A 
transition from an activated regime to an entropic one is observed at temperature Tg /2 in agreement 
with recent work on this model. Ergodicity breaking occurs for T < Tg/2 in the thermodynamic 
limit, if K/N 0. In this temperature range, the model exhibits a non trivial fluctuation-dissipation 
relation leading for K <^ N to a, single effective temperature equal to Tg/2. These results give new 
insights on the relevance and limitations of the picture proposed by simple trap models. 



I. INTRODUCTION 

An important step towards the understanding of glassy 
dynamics has been made when it was recognized that 
some generic properties of configuration space -or phase 
space- could be responsible for the dramatic slowing 
down of the dynamics ^^J. In particular, the ge- 
ometric structure of phase space leads schematically to 
two different kinds of dynamics: an 'activated' dynam- 
ics in which the system is trapped in local minima by 
significant energy barriers, and an 'entropic' dynamics 
which results from a decreasing number of downwards 
directions when visiting saddles in configuration space 
IE 1^ Qi ^1 ' In this latter case, the system spends 
most of its time wandering in search of these rare paths 
which would allow it to decrease its energy. 

A popular and qualitative description of these glassy 
behaviors has been proposed in the past decade in terms 
of trap models, in which a very simplified phase space 
dynamics takes place. In these models, any state can 
be reached from any other through a single transition, 
disregarding any non trivial structure related to the finite 
dimensionality of real space. Such models actually focus 
on the distribution of low energy states, often assumed 
to be exponential, which can be justified on the basis of 
extreme statistics [T^ . 

Depending on the specific choice of the transition rates, 
one can build an activate d dy namics -as in Bouchaud's 
trap model (BTM) [ijjIlJjO" or an entropic one -as in 
the Barrat and Mezard model (BMM) [Ulil. Consider- 
ing a finite size BMM, or introducing by hand a threshold 
level, one can observe a crossover from an entropic to an 
activated regime |l7j . Intuitively, such a crossover means 
that the system is no longer able to find downwards di- 
rections since it has reached the bottom of the 'valley'. 
Further evolution can only proceed by crossing energy 
barriers. 

In spite of the conceptual interest of these models, 
it seems rather difficult to find microscopic models (i.e. 
models in which microscopic degrees of freedom are ex- 



plicitly described) where a reasonably clear mapping to 
such trap models can be proposed. This situation is par- 
ticularly striking since the physical interpretation of trap 
models looks quite clear, but arguments usually fail to 
go beyond a qualitative level. 

The first explicit (and mathematically rigorous) map- 
ping was proposed between the finite size Random 
Energy Model [I^ and the BTM. Trap mechanism has 
also been shown to be a tangible description of super- 
cooled liquids slowing down when considering the distri- 
bution of the energy associated to the inherent structures 
[2y|. On the other hand, it has been proposed recently 
|2l| to use a modified version of the Number Partition- 
ing Problem (NPP), mapped onto a fully connected spin 
model with a one-spin-flip dynamics, to illustrate how 
an activated behavior typical of the BTM arises from a 
microscopic dynamics. 

In the present paper, we discuss the influence of the 
choice of the dynamics on the behavior of the NPP. We 
show in particular that varying the number of spins that 
can be flipped simultaneously allows to recover most of 
phenomenology of glass theory, namely transitions be- 
tween entropic and activated behavior, non-linear as well 
as linear (with non-trivial slope) fluctuation-dissipation 
relation (FDR), and ergodicity breaking. Conversely, 
these microscopic realizations allow to shed some new 
light on the interpretation, as well as limitations, of sim- 
ple phase space models like the BTM and the BMM. 

The paper is organized as follows: Sect. II introduces 
the NPP model, and Sect. Ill describes the basic map- 
pings onto usual trap models for some speciflc dynamical 
rules. In Sect. IV, we introduce more general dynamical 
rules, and study the behavior of the model, emphasizing 
the entropic-to-activated transitions as well as relations 
to trap models. In Sect. V, the FDR is studied and shown 
to be linear in a particular limit, with a non trivial slope. 
Finally, we discuss in Sect. VI the interpretation of this 
linear FDR, as well as the influence of the energy density 
on the transition between entropic and activated regimes. 
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II. OPTIMIZATION AND SPIN MODEL 

The optimization problem of the unconstrained NPP, 
described by a cost function, can be mapped onto a f ully 
connected spin model with disordered Hamiltonian [23, 
01 ■ By suitably choosing the cost function and the 
dynamics, this spin model can be given a glassy behavior 
which resembles closely that of the BTM 21J. Given 
a set of N random real numbers 01,02, ...aN uniformly 
drawn from the interval [0, 1], the original NPP consists 
of finding the optimum configuration {si}°^-^ ^, where 
Si — ±1 are Ising spins, which minimizes the following 
cost function: 
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This is equivalent to finding a partition of the set {oi} 
into two subsets Si and S2 such that the sums of the 
fli's within each subset are as close as possible. In terms 
of spin systems, this problem corresponds without loss 
of generality to an Hamiltonian E'lviattis = E'^ describing 
an antiferromagnetic Mattis-like spin-glass: 
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From a thermodynamic point of view, Mertens j22, |23| 
has shown that the ground state of this Hamiltonian was 

{Em J = y1 



so that no glass transition at finite 

temperature is expected in the thermodynamic limit. In- 
terestingly, from such an Hamiltonian, one can derive a 
new cost function, i.e. a new energy, that has an extensive 
ground state: 
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where Tg fixes the energy scale; from now on the ground 
state scales like (Eo) ^ —NTg\n2. In this paper, we 
consider the study of a system defined by such an Hamil- 
tonian. 

In this system, energies E < A; In TV (where k is some 
positive constant) are independent random variables pH 
12^ that are distributed according to (assuming N ^ 1): 



p{E) = Af exp ( PgE 
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with Af = 2Pg/V2na'^N, and (3g = Tg^^. The essential 
property of this distribution p{E) is that it has an expo- 
nential tail for E' — > — 00. Such a tail is usually the key 
ingredient to obtain a glass transition at finite tempera- 
ture. 



Using Derrida's microcanonical argument for the ran- 
dom energy model (REM) _19J, a thermodynamic transi- 
tion is expected at temperature Tg , below which the sys- 
tem is frozen in a limited number of states surrounding 
the ground state, so that the entropy (density) vanishes. 

The glass transition in the present model resembles the 
standard REM transition, the only difference being that 
the former is first order whereas the latter is second order 
0,1^. In particular, an important property shared by 
both models is that for low energy states, magnetization 
and energy become decorrelated. 

From an optimization point of view, many interesting 
questions are inherent to the NP-hard nature of the NPP 
lii,!!^. Thus, it is interesting to study, when T < Tg, 
how the system (j^l approaches the ground state depend- 
ing on the local dynamical laws, given the prescription 
of detailed balance. Since the NPP belongs to the class 
of NP-complete problems, time needed for any algorithm 
to get the perfect partition is exponential in the system 
size N. The most naive algorithm which consists of an 
exhaustive enumeration of all the partitions is then as 
efficient as any elaborated algorithm when N becomes 
large [2^ . Subsequently, even though the dynamics stud- 
ied in this paper seems a priori to be inappropriate to the 
optimization problem, all the underlying mechanisms re- 
sponsible for out-of-equilibrium processes, especially ag- 
ing phenomena, come from the NP-hard nature of the 
problem. 

So below Tg any dynamical rules (satisfying detailed 
balance) lead to an aging regime before reaching the 
ground state. In the following, we use K-spin-flip 
Metropolis rules {I < K < N) defined as follows. At 
each Monte-Carlo time step, a new configuration {s'^} 
that differs from the current one {si} by at most K spins 
(and at least one) is chosen randomly. This new config- 
uration is then accepted with a probability equal to the 
Metropolis acceptance rates at temperature T: 



p-{E'-E)/T 
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Monte-Carlo time steps are separated by a physical time 
interval r,„c — K/N in the natural time units of the 
system. This ensures that each spin keeps a probability 
of the order of 1 to be chosen within a unit time interval, 
even in the thermodynamic limit. 

We show in the following that this model leads to a 
rich dynamic phase diagram, the control parameter being 
the fraction K/N of spins allowed to flip simultaneously. 
This phase diagram can be discussed in light of both the 
BMM and BTM studies at finite temperature. 



III. SIMPLE REALIZATIONS OF TRAP 
MODELS 

In this section, we study two different dynamical rules: 
a single-spin-flip dynamics {K = 1) and a global dynam- 
ics involving full rearrangement [K = N). Interestingly, 
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FIG. 1: Schematic representation of the phase space structure 
of the NPP model with single-spin-flip dynamics. The hori- 
zon level separates surface states and low energy states. An 
observable is smooth if it varies slowly between neighboring 
states. 



these two limiting cases appear to be microscopic real- 
izations of trap models, the former with an activated be- 
havior and the latter with an cntropic one. 




FIG. 2: Aging behavior of the correlation function C{tw,tw + 
t) in the NPP with /C = iV = 50, for T = 0.75Tg and 
different t-u); data rescale as a function of t/t^. The full 
line is the analytical prediction Eq. IIUII for t ^ t^. In- 
set: correlation for different temperatures; from left to right: 
T/Tg = 0.75,0.6,0.5,0.4,0 = 10^). The dashed line is 
proportional to . 



A. Single-spin-flip dynamics: activated traps 

It has been shown recently that a single spin flip 
dynamics naturally leads to an activated trap behavior a 
la Bouchaud, in which the aging phenomenon comes from 
the divergence of the average trapping time. Given the 
density of state l@J with an exponential tail, two dynami- 
cal ingredients are responsible for such a behavior: on the 
one hand, the existence of an horizon level below which 
the system has no choice but to reemerge above it so as to 
continue its evolution; on the other hand, instantaneous 
jumps into a randomly chosen new trap after reemerging 
at the horizon level, associated with a full decorrelation. 

The former appears naturally since when the energy is 
lower than: 

Eh=Tg\n{a,r.^)K-Tg\nN (6) 

with flmin = min(ai, ajsi) = 0{1/N), a single spin flip 
necessarily leads to a state whose energy is greater than 
Eh- The latter is due to the combination of two prop- 
erties. First, low energy states are totally uncorrelated; 
second, the travel time around Eh between low energy 
states becomes negligible as time increases, that is as 
lower and lower energy states are visited. 

Interestingly, the need for reorganization around high 
energy levels in order to go from one low energy state 
to another is responsible for an equilibrium-like linear 
FDR with slope 1/T for smooth observables, i.e. observ- 
ables like the magnetization whose relative variation is 
of the order of \/N after one spin flip -see Fig. ^ for a 
schematic view. This law is observed even in the aging 
regime, T being the temperature of the thermal bath. 



This comes from the fact that the evolution of smooth 
observables is dominated by the sojourns among high en- 
ergy states, where they can (almost) equilibrate. Once in 
a deep state, these observables become frozen, but their 
typical value is indeed that given by Gibbs distribution, 
since it is determined by the high energy states visited 
just before falling into the trap. 



B. A'^-spin-flip dynamics and BMM behavior 

Let us consider now a global dynamics (i.e. K = N) 
such that all spins are flipped randomly (and simulta- 
neously) at each step in order to find a new configura- 
tion. The transition is then accepted or rejected accord- 
ing to the Metropolis rates given in Eq. (|SJ). As a result, 
any configuration is a priori accessible from any other 
(apart from the energetic constraint), which means that 
the horizon level disappears, and the new configuration is 
in general completely decorrelated from the old one. As 
moreover energies are distributed exponentially, one can 
expect this model to be a microscopic realization of the 
BMM. In the following, we propose more quantitative ar- 
guments as well as numerical simulations to support this 
statement. 

In all the numerical simulations, we have dealt with 
the autocorrelation function C {tw ,tw + t) between time 
and + t defined by the average over the thermal 
histories of the history-dependent correlation function 

C Single {^w : t) . 

+ t)) (7) 
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FIG. 3: Short time behavior of the correlation function 
C{tw,tw+t) and onset of a singularity with exponent (1—^) //i 
(fuU hnes) for § < fi< 1 {ii = T/Tg), as given by Eq. |[TH 



with 



1 ^ 

C Single ij^w i H" ^) = "J^ ^ ^ Si {tw^Si [tyj + ^ 



(8) 



This choice of correlation is usual in spin models, al- 
though other choices like the autocorrelation of the mag- 
netization would be possible. Actually, C{tu,,tu]+t) is the 
autocorrelation of the observable X^iCi^i, where = ±1 
are quenched random variables. We show in the appendix 
that the specific choice of the observable does not influ- 
ence the main properties of the model, as long as the 
observable is smooth. 

In the case of a full redistribution of the spins, this 
autocorrelation reduces to the hopping correlation func- 
tion C^{tw,tw + t), defined by the following history- 
dependent function: 



+ t) = 



if 



otherwise 

(9) 

which precisely leads to the same correlation as in trap 
models. The aging regime is characterized by the fact 
that the correlation function C{tw,tw+t) becomes a func- 
tion C{t/tw) of the ratio i/tu, only. For trap models with 
exponential energy distributions, the asymptotic behav- 
ior of the function C(t/tu,) for t <^ tw and for t ^ t^, is 
characteristic of the nature (entropic or activated) of the 
dynamics Tt'I. In the NPP, numerical data show that 
(full) aging is observed for any temperature T < Tg as 
expected (see Fig.|21for T = 0.75Tg). 

In the long time limit {t 3> i^i), the asymptotic be- 
havior characteristic of the BMM, with a temperature 
independent exponent, is recovered (Fig. [SJ: 



Actually, one can be more specific and compute the 
exact asymptotic expression of the correlation function 
in the case of Metropolis rates. One finds C{tw,tw +i) ~ 
{1 — fj,)tyj/t, where fi — T/Tg is the reduced temperature. 
This prediction fits well the numerical data, as shown on 
Fig. El 

In the short time limit (i <C t^,), an asymptotic anal- 
ysis of the correlation function in the BMM shows the 
onset of a singularity for temperatures in the range 
Tg/2<T<Tg: 



1 C(tjjj , til] + 1) 



(i) 



i < /i < 1 



(11) 



(10) 



More precisely, the singularity concerns the scaling func- 
tion C{u) in the limit u — *■ 0. This singularity clearly 
appears also in the NPP with the A^-spin flip dynamics, 
as shown on Fig. |3| This is the signature of an entropic- 
to-activated transition at temperature Tg/2 Dis- 
crepancies at very short time come from an exploration 
of states that are not exponentially distributed due to 
finite size effects, but also from finite time effects since 
correlation functions are calculated analytically in the 
limit of asymptotically large times. 

Different kinds of transitions between entropic and ac- 
tivated dynamics have been found in the context of the 
BMM, or of modified versions of this model One, al- 
ready mentioned in the introduction, is a crossover from 
an entropic to an activated regime as the system ages. 
The heat bath temperature is kept constant in this pro- 
cess, and a characteristic time scale is associated to the 
crossover. A second type of entropic to activated transi- 
tion can be found also when varying temperature. In the 
BMM with exponential energy density, such a transition 
appears when temperature crosses the value Tg/2] below 
Tg/2^ the dynamics is essentially dominated by entropic 
effects, while above this value (with still T < Tg to re- 
main in the aging regime) activated effects come at play. 
This is seen in particular from the short-time behavior 
(i <C iu)) of the aging correlation function which becomes 
singular above Tg/2^ with an exponent (1 — /i)//i remi- 
niscent of (although different from) the exponent 1 — /i 
of the BTM This transition is also present in the 

NPP with K = N, as seen from Eq. Hll|l . confirming the 
mapping between both models. 



IV. INTERMEDIATE DYNAMICS: 1< if < iV 

We have seen in the previous section that the limit- 
ing dynamical rules {K — 1 and K — N) correspond 
to the two simple kinds of trap models, namely BTM 
(activated) and BMM (entropic). Since the ratio K/N 
governing the dynamical rules can be varied (almost) con- 
tinuously from to 1 , a crossover between both kinds of 
behaviour should be found. Yet, this crossover is rather 
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non trivial, as activated and entropic dynamics are qual- 
itatively different. In particular, one can expect the hori- 
zon level to play a major role in this change of dynamics. 
The way the observables decorrelate may also lead to 
significant differences with respect to the standard trap 
picture. 

A. Differences with the previous rules 

1. Slow decorrelation of smooth observables 

In simple trap models, one usually assumes that the 
value of the observable after a transition (a 'jump') is 
completely decorrelated from the value it had before the 
jump. This simple assumption, which might seem unreal- 
istic at first sight, is indeed fulfilled by the single-spin- flip 
and the N-spin-flip dynamics, but for different reasons. 
li K = N, it is clear that at each step, the new configu- 
ration is independent of the old one, so that observables 
immediately decorrelate. For the case K = 1, the mecha- 
nism appears to be more subtle: each time a new configu- 
ration is chosen, smooth observables typically decorrelate 
by a factor (1 — 1/N). Dynamics is dominated by low en- 
ergy states (or traps), which are below the horizon level. 
Once in such a trap, a single spin flip leads to a high 
energy state, and many subsequent flips are necessary in 
order to find a new trap. Yet, the typical time spent wan- 
dering among these high energy states remains negligible 
compared to the time spent within traps. So in terms 
of the effective dynamics between traps, the observables 
indeed fully decorrelate at each jump. 

So what happens for 1 ^ if ^ A^? In this case, two 
subsequent configurations are still highly correlated, and 
it is not clear either whether some relevant coarse-grained 
description could lead to an effective decorrelation. One 
thus expects to observe some non trivial behaviors which 
may differ significantly from that of usual trap models. 

2. Influence of the horizon level 

The single-spin-flip dynamics studied previously has 
emphasized the fundamental role played by the horizon 
level. On the other hand, this threshold completely dis- 
appears for K — N. In the intermediate case 1 ^ K 
N, one can still define an horizon level, but this level is 
expected to drift towards lower energies as K increases. 
Using the same argument as above to determine the hori- 
zon -see Eq. one gets the following threshold energy: 

= -^KTg In N {K < N) (12) 

Below this level, the evolution is always activated: an 
energy barrier at least equal to {Ef^ — E) has to be over- 
come when starting from an energy E < E^ . 

On the contrary, as long as the system visits states with 
energy E well above E^ , the influence of the threshold 
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FIG. 4: Schematic view of tlie different dynamical regimes 
appearing in the NPP depending on T and on the typical 
energy E visited by the system (large values of \E\ correspond 
to deep energies and long times). 

should not be felt. So one can guess that two differ- 
ent dynamical behaviors for temperatures above and be- 
low Tg/2 should still exist, as found also in the modified 
version of the BMM including a threshold ^| . As a re- 
sult, one expects to find schematically the three following 
regimes (see also Fig. ^: 

1. E > £'f , T < Tg/2. This case resembles the en- 
tropic regime of the BMM, the only difference be- 
ing that the magnetization decorrelates slowly, typ- 
ically by a factor (l—K/N) after each jump, at vari- 
ance with the full decorrelation usually assumed in 
the BMM. 

2. E > Ejf, T > Tg/2. As in case (1), the dynamics 
is similar to that of the BMM; in this temperature 
range, activated effects become important, and ob- 
servables also decorrelate slowly. 

3. E < E^ . Once the horizon level is reached, the sys- 
tem has no choice but to reemerge above it, which 
leads to an activated dynamics similar to that of 
the BTM. No specific role is played by the temper- 
ature Tg/2 below the horizon, this activated regime 
is qualitatively the same in the whole range T < Tg. 

Thus, the NPP with a finite number of spin fiips 
(i.e. K <C N) exhibits a crossover from an entropic regime 
to a activated one as time elapses, for any given temper- 
ature T < Tg. Although it would have been interesting 
to investigate the properties of the model within (and 
beyond) this crossover regime, we do not study them in 
the present paper, since the crossover time tx is exponen- 
tial in K (more specifically, ~ exp{-El^ /Tg) ~ N'^) 
which leads to time-consuming numerical simulations. 
We thus focus on cases (1) and (2), corresponding to en- 
ergies well above the horizon level. Subsequently, all the 
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results reported below correspond to times much smaller 
than the crossover time t y . 



B. Entropic versus activated dynamics 

1. Qualitative approach 

Before giving quantitative arguments about correlation 
functions, let us propose a more intuitive understanding 
of the difference between entropic and activated dynam- 
ics. To this aim, it is of interest to plot the energy as a 
function of time for a single thermal history. This is done 
on Fig. [SI for three different temperatures, T/Tg = 0.35, 
0.5 and 0.65. For T < Tg/2, the energy decreases es- 
sentially in a monotonic way, and the evolution is close 
to the zero temperature one: the dynamics remains en- 
tropic. On the contrary for T > Tg/2, the system comes 
back many times to high energy levels, as if it had to 
reemerge from a deep trap: activated events dominate 
the evolution, which is then rather similar to that of the 
activated trap model. 

Note that qualitatively similar trajectories in energy 
space can be found in the case K = N discussed in the 
previous section as well as in the BMM. In this context, 
it has been proposed 0] to characterize the type of dy- 
namics by computing the average energy {E')e reached 
in a transition between two different microscopic states, 
starting from a given energy E: 



E') 



r dE' E' WiE 

{E')e = ^ 

/° dE'W{E^E') 



(13) 



Using the Metropolis rules and assuming that the vis- 
ited energies are still well above the horizon level, this 
quantity reads in the large \E\ limit: 



{E')E = E+^J^Tg 



(14) 



Below Tg/2, the energy is lowered on average at each step, 
leading to an irreversible drift towards low energies. 

On the contrary, above Tg/2 the energy is raised on 
average, i.e. {E')e > E, as long as \E\ is large enough 
{E < 0). So the energy variable performs a random walk 
in energy space, with a bias towards high energies, and 
-roughly speaking- a reflecting boundary condition in 
E = 0. If time was counted in number of jumps, the 
random walk would then reach a steady state: energies 
tend to remain close to the boundary E = 0, as seen on 
Fig. El However, the larger \E\, the larger the sojourn 
time, so that large fluctuations away from the boundary 
(i.e. at large \E\) dominate the real time dynamics. So 
the probability to have energy E at time t does not reach 
any steady state, and drifts continuously towards low en- 
ergies. This competition between the bias towards high 
energy and the large trapping time of low energy leads 
to the onset of a singularity in the correlation function 
as discussed above. 
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FIG. 5: Energy as a function a time for a single thermal 
history at three different temperatures: from left to right, 
T/Tg = 0.35, 0.5 and 0.65 {N = 200 and K = 5). Returns 
to high energy levels appear only for T > Tg/2. At this 
temperature, the dynamics is rather similar to the BTM, with 
in particular an apparent reversibility when plotted on a linear 
time scale I21I. 



Since a full analytical solution of the correlation func- 
tion in the present model with K <^ N seems difficult to 
reach, simple scaling arguments can be helpful in order 
to interpret numerical simulations. Data shown on Fig. [3 
suggest that the type of scaling argument (or in other 
words, the relevant approximations) may be different for 
entropic and activated dynamics. In the entropic range 
T < Tg/2, the instantaneous energy E{t) can be decom- 
posed into an average value (with deterministic evolu- 
tion) and a fluctuating term with zero mean and a finite 
amplitude. Even though fluctuations are not necessarily 
small, they do not dominate the dynamics and may be 
considered as a perturbation over the average determin- 
istic evolution. So it may be reasonable to think that 
a zeroth order approximation which would neglect fluc- 
tuations could yield some relevant results, in particular 
concerning the scaling behavior. 

On the contrary, the dynamics for T > Tg/2 appears 
to be dominated by activated events during which the 
system visits high energy states. Fluctuations are now 
driving the evolution, and cannot be considered anymore 
as a perturbation which could be ignored in a first step. 
Indeed, since the system goes back frequently to super- 
ficial states at energy E k, Q, the amplitude of the fluc- 
tuations (with respect to the average energy at time t) 
diverges with time. So scaling arguments involving only 
average values cannot be used anymore. 



2. Entropic temperature range T < Tg/2 

Let us first consider the case of zero temperature and 
estimate the aging law for magnetization through a sim- 
ple scaling argument. Given that the magnetization 
decorrelates typically by a factor {1 — K/N) at each tran- 
sition, one can compute the correlation Cr after i? jumps. 
Assuming K <^ N, we have in the large R limit: 



-RK/N 



(15) 



with the prescription Cji=o = 1. From Eq. (|14|l . we know 
that after each jump, the energy decreases of Tg on aver- 
age, so that after R jumps the energy difference between 




FIG. 6: Correlation function C{tw,tw + t) up to very long 
time t at zero temperature and for small values of the ratio 
K/N . Data were obtained from an efficient (event-driven) 
Monte-Carlo algorithm; the thermodynamic limit A'' ^ cx) is 
taken, so that the horizon level = —KTglnN cannot be 
reached. 
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FIG. 7: Direct Monte-Carlo simulations of the correlation 
function in the range T < Tg/2, for small values of K/N. For 
T > 0, iV = 200 and U = 10*; for T = 0, = 1000 and 
tw ~ 500. It has been checked that typical energies remain 
well above the horizon level. Lines are predictions given by 
Eq. I22II : no fitting parameter is used. 



times and + 1 is given by: 



(16) 



From an energetic point of view, the NPP far above Ej^ 
is expected to be equivalent to the BMM. Subsequently, 
at an energy E the corresponding trapping time is given 
by: 



TE = T„ 



dE'p{E') 



(17) 



where p{E') is given by Eq. Q and Tmc = K/N -see 
Sect. ^ At low energy, this trapping time reduces to: 



r{E) 



(18) 



where A/" is a factor coming from the distribution p{E). 
Since time -I- i is of the order of the typical trapping 
time of the state currently visited, Eq. H18|l leads to the 
following relation between t^, -I- i and E: 



N Af 

This combined with Eqs. ((T3|l and gives: 



-K/N 



(19) 



(20) 



for the aging correlation function at zero temperature. 
Although rather naive, this simple scaling argument is 
very well confirmed by long time simulations performed 
with an efficient event-driven algorithm, as seen on Fig.El 



The agreement with direct Monte-Carlo simulations of 
NPP up to accessible times is also very satisfactory (see 
Fig. lower curve). Note that a modified version of the 
BMM which includes the property of slow decorrelation 
of the observable precisely leads to the same behavior as 
Eq. H20() at zero temperature '2^. 

For non zero but low enough temperature (T < Tg/2), 
so that activated processes do not dominate the dynam- 
ics, one can use exactly the same argument, simply mod- 
ifying the time dependence of the average energy accord- 
ing to Eq. ifTH) : 



E{t^, + t) - E{t^ 
This gives: 



-^RTg (21) 



tin 



-)] K/N 



(22) 



with 1] = (1 — /i)/(l — 2/i). Here again, this simple es- 
timation describes rather well the numerical simulations 
(Fig.Q. 

So one sees that the law of decorrelation of the ob- 
servable between states, which is encoded in the ratio 
K/N, has a dramatic impact onto the corresponding ag- 
ing laws. In particular, the specific behavior of the cor- 
relation function given by Eqs. ^ and ^ has im- 
portant consequences regarding the thermodynamic limit 
N 00. According to the dependence of K on N, 
the correlation function is able or not to decay at large 
times. Indeed, if X w aN with some positive constant 
a, C(t^,tw + 1) converges in the large TV limit to a well 
defined scaling function which decays to for t ^ 00 -see 
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FIG. 8; Correlation function in the regime Tg/2 < T < Tg 
for the NPP with 1 < A" < A^. Straight lines have slope fi = 
T/Tg on log- log scales, so as to compare with the activated 
behavior of the BTM. 



Eq. (Unj). On the contrary, if K/N ^ for TV ^ oo (say 
K is fixed), the system becomes unable to decorrelate in 
the thermodynamic limit, and C{tw, tw+t) remains equal 
to 1. Defining a as the limit when iV ^ oo of the ratio 
K/N, and taking it as a control parameter, one sees that 
a transition towards a state where ergodicity is broken 
occurs at a = 0. 



3. Activated temperature range T > Tg/2 

As already mentioned in the introduction of this sec- 
tion, the dynamics in the temperature range Tg/2 < 
T < Tg is qualitatively different from that in the range 
T < Tg/2: jumps to high energies, inducing large energy 
fluctuations, play an important role in the evolution of 
the system. In this case, a scaling argument based only 
on the deterministic evolution of average values is not 
expected to be relevant. 

Let us consider first the dynamics of the energy. Since 
for K 1 the horizon level Ej^ is low, a significant 
energy range above Ef^ exists where energy states are 
all mutually accessible, independent and exponentially 
distributed. Given the dynamical rules Eq. (0, the evo- 
lution of the energy is precisely that of the BMM; for 
instance, one should find the same dynamic probability 
distribution P{E,t) -i.e. the probability to have energy 
E at time t. Besides, once the threshold level is reached, 
a fully activated dynamics typical of the BTM should be 
recovered. So as far as the evolution of the energy is con- 
cerned, the situation is very similar to the BMM with 
threshold studied in [T^ . 

Turning to the correlation function of smooth observ- 
ables, the situation is a bit more subtle. Numerical data 
from the NPP model are shown on Fig. |S1 At variance 
with the usual results of the BMM, the long time tail 




FIG. 9: Correlation function in the BMM with a decorre- 
lation factor (1 — x); here, x — 0.1. From top to bottom: 
T/Tg = 0.52, 0.55, 0.6, 0.7 and 0.8. The long time tail be- 
comes temperature dependent, with an exponent close to the 
value — /i found in the BTM, but deviations increase when 
T Tg/2. Inset: short time behavior, with the typical expo- 
nent (1 - fi)/fi of the BMM (dashed). 



(t ^ tw) of the correlation function seems to behave as a 
power law with a temperature dependent exponent. Yet, 
one must admit that these numerical data do not allow 
to draw a definitive conclusion on this point. 

Since the magnetization decorrelates only by a factor 
(1 — K/N) (with K ^ N) at each transition, the evolu- 
tion of the correlation should be that of the BMM with 
slow decorrelation. Still, this model has not been studied 
in details in the literature, and one could wish to know 
whether the behavior of the correlation function resem- 
bles that of the hopping correlation function in one of the 
usual trap models. 

To this aim, we have simulated directly a modified ver- 
sion of the BMM in which the observable is decorrelated 
only by a factor (1 — x) at each transition, assuming 
X <C 1. Numerical data are shown on Fig. |51 Interest- 
ingly, the behavior of the correlation function is numeri- 
cally found to be reminiscent of both the BTM and the 
BMM. The short time behavior (t <C tyj) is very similar to 
that of the BMM: a singularity 1 — C ~ {t/t^)'* appears, 
with an exponent 7 very close to the value {l — fi)/fi found 
for the hopping correlation function of the BMM. On the 
contrary, the long time tail becomes temperature depen- 
dent, at variance with the usual BMM behavior given by 
Eq. (|10|) . and in agreement with numerical results found 
in the NPP -see Fig. |S1 The corresponding exponent is 
close to the exponent — /i typical of activated regimes, 
but significant discrepancies appear for T close to Tg/2. 
These discrepancies were somehow expected from a conti- 
nuity argument, since below Tg/2 the correlation function 
decays with a very small exponent —rjK/N . 



V. FLUCTUATION-DISSIPATION RELATIONS 

Generalizing the well-known equilibrium theorems, 
fluctuation-dissipation relations (FDR) have proven to be 
a very robust tool to study the out-of-equilibrium regime 
of glassy models ■ Given an observable O, in most 
cases |23i|23|, the two points correlation function between 
time ti and t2 > ^i, 



(23) 



and the response to a perturbation h conjugate to O, 



are related through a FDR: 

7^(^l,<2) = 



Sh{ti) 



h=0 



dC{hM) 



(24) 



(25) 



At equilibrium, Te//(ti,i2) is given by the thermal bath 
and all the system properties are time translational in- 
variant. Far from equilibrium, these FDR can be used 
to define a meaningful effective temperature [2^ |2^ , 
as in the context of mean-field (or fully connected) spin- 
glass models • Since we are indeed dealing here with a 
fully connected spin model, it is then natural to study the 
FDR. An important question to address is the temporal 
independence of Teff{ti,t2) since only when Teff(ti,t2) 
does not depend on times can a unique and well-defined 
effective temperature be introduced. In this case, intro- 
ducing the integrated response, 



dt'n(t',t2) 



(26) 



the FDR is said to be linear since the i2-parametric plot 
x(ti, t2) vs. C(ti, t2) is a straight line whose slope is given 
by -l/Teff. 



A. FDR in the aging regime for K — N 

The temperature Tg/2 separates the two different 
classes of dynamics encountered in this model. It is then 
of interest to compare the FDR in these two regimes. Nu- 
merical data concerning the FDR for the global dynamics 
{K = N) are shown on Fig.jTUl in the case T < Tg/2. 

As the NPP model with a global dynamics can be 
mapped onto the BMM, one expects to recover the exact 
result found in the BMM at zero temperature [s^ -see 
also |33|- which reads: 



t) 



2 dC{U 



t) 



Tn 



dt 



(27) 



R{tw+t, tu;) being the response associated to the autocor- 
relation function C{tw, tm + t). In terms of the integrated 
response, this FDR can be reformulated as: 

x{tw,t^+t) = ^[l-C{t^,t^+tf] (28) 




C(t ,t +t) 

FIG. 10: Fluctuation-dissipation relations in the NPP for 
K = N — 50 for different temperatures (t^ = 10'*); x is 
measured in units of Tg. Numerical data at T = Tg/2 remain 
very close to the zero temperature analytical relation of the 
BMM (full line). For T > Tg/2, the initial slope is given by 
l/T (dashed line). 



taking into account the explicit expression of the corre- 
lation [s^l- Numerical data shown on Fig.^jare in good 
agreement with this analytical prediction. Interestingly, 
this zero temperature solution seems to be also valid 
for any temperature in the 'entropic range' T < Tg/2 
fFig. llO|l . As a result, below Tg/2, the FDR remains non 
linear. 

From a technical viewpoint, the autocorrelation (9) is 
associated to the observable where {^i} is a 

set of quenched independent random variables that can 
take the values ±1. It is understood that all measured 
quantities are averaged over the realizations of {^i}. So 
the integrated response xitwitw + t) is numerically ob- 
tained in the NPP model by computing 



1 ^ 

xit^„t^+t) = —Y^,s^{t^+t) t>0 (29) 



i=i 



where /i is a small external field that is switched on at 
time t^. This field h is coupled to the spins via a linear 
coupling term Vh in the energy, involving also ^i: 



N 



(30) 



Above Tg/2, the FDR significantly depends on T but 
seems to have unchanged behavior in the two following 
limits: t <^ and t ^ tw (Fig. ITUI) . In the former, 
the slope of the curve C vs. R is T^^ (as in equilib- 
rium) whereas in the latter, the slope apparently goes 
to 0, which might be interpreted as an infinite effective 
temperature. Yet, such a temperature must be taken 
with care since the definition of effective temperatures 
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FIG. 11: FDR for T > Tg/2 in the aging regime of the NPP, 
for 1 <^ K <^ N: a linear relation with slope is observed. 
Simulation parameters: A'' = 1000 and tw ~ 10^ for T — 
Q.GTg- N = 100 and i„ = 2 x lO'^ for T = 0.7Tg. 

from the local slope of non-linear FDR remains an un- 
clear procedure. 

Apart from this question of effective temperature, it is 
quite interesting to see that one can discriminate between 
an entropic regime and an activated one in the BMM at 
finite temperature, from the initial slope of the FDR. 
The entropic regime gives a T-independent slope corre- 
sponding to the temperature Tg/2 that separates the two 
regimes, whereas the activated regime gives a slope that 
corresponds to the thermal bath temperature T. 



B. Out-of-equilibrium FDR for < 

Let us consider now the intermediate dynamics 1 ^ 
K <^ N. In the case T > Tg/2, one recovers when 
K/N — > the same behavior as for X = 1 namely 
a linear FDR with a slope equal to 1/T, where T is the 
heat bath temperature. Thus there is no violation of 
fluctuation-dissipation theorem, as seen on Fig. ^2 The 
mechanisms at play are essentially the same as for the 
activated dynamics found when K = 1 -see 113" so 
that this case need not necessarily be discussed in de- 
tails. Note that although simple trap models usually do 
not display linear FDR [s^ [33, 13 1 such linear relations 
with slope 1/T as in equilibrium can indeed be found 
in trap models (at least for some particular observables) 
once a spatial structure is taken into account 35] . 

On the contrary, the behavior of the system for T < 
Tg/2 is much more surprising. Indeed, as K/N 0, 
numerical data on the zero temperature FDR converge 
to a linear relation with slope 2/Tg. Results for differ- 
ent values of K/N and different temperatures below Tg/2 
are gathered in Fig. E| As mentioned above, the limit 
of linear FDR allows to define an effective temperature, 



FIG. 12: FDR in the entropic temperature range T < Tg/2 
in the NPP with 1 <^ K <^ N . Data were obtained by direct 
Monte-Carlo simulations (MC) or by the event-driven algo- 
rithm (ED); tyj = 10^ for MC and 10^° for ED. Data converge 
to a linear FDR in the limit K/N ~* (and — » oo), with a 
slope 2/Tg. Thus in the regime, the effective temperature is 
equal to Tg/2. 



which would be equal here to Tg/2 and is thus tempera- 
ture independent. So it appears that this transition tem- 
perature between activated and entropic dynamics again 
plays an important role in the description of the low tem- 
perature phase. It is also worth noticing that the initial 
slope (i.e. for C ~ 1) of the non-linear FDR found for 
the global dynamics K — N are the same as the slopes 
of the hnear FDR in the case 1 < X < TV. 

Here again, in this regime T < Tg/2, the NPP model 
presents strong similarities with the BMM modified to 
include slow decorrelation, as discussed in Sect. IIV B 21 
Indeed, it can be shown that this particular version of the 
BMM, with a vanishing decorrelation of the observable 
at each transition, also leads to a linear FDR with a slope 

Tg/2 m 



C. Interpretation of the linear FDR 

Interestingly, this linear FDR with slope 2/Tg can be 
derived analytically in the whole regime T < Tg/2 for 
smooth observables. The corresponding calculations are 
reported in Appendix. In this section, we simply try 
to give a simplified picture of the physical mechanisms 
leading to this non-trivial effective temperature by con- 
sidering the case of magnetization. 

As seen in the case K — 1 ,2lj, a key to the linearity 
of the FDR is the fact that the magnetization induced 
by the field h is not influenced by the state in which the 
system was at time tw, when the field was switched on. In 
other words, the contribution to the magnetization of this 
initial state should be negligible. This is indeed the case 
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here, since the contribution is expected to be of the order 
of Kh, whereas the total magnetization should scale as 
Nh. So in the limit K/N 0, the above contribution 
vanish. Notice that this property does not hold for the 
BMM (A' — N) since the initial state contribution is 
expected to be of the order of Nh. Subsequently, this 
case is not expected to give a linear FDT in agreement 
with the established results ^32,] (see Eq. (|27|l ). 

Once this contribution is neglected, the linearity of the 
FDR can be given through a rather simple physical inter- 
pretation. It simply means that the relaxation towards 
the non-zero magnetization induced by the field h be- 
haves in the same way as the relaxation towards zero 
starting from an arbitrary magnetization induced by the 
spontaneous fluctuations. So the slope of the FDR is de- 
termined by the asymptotic value of the magnetization, 
a long time after the field was applied. For K = 1, the 
visits to superficial states induce asymptotically a con- 
stant magnetization M*, since the a priori distribution: 



p{M) 



(31) 



is weighted by the Boltzmann factor e^^^l^ \ the average 
value of the resulting distribution is then M* = hN/T. 
On the other hand, the equal-time correlation (in the 
absence of field) C{t,t) = ( M'^) /N is equal to 1, hence 
the slope 1/T of the FDR [2J. One can expect this 
argument to be valid also in the regime T > Tg/2 for 
1 ^ A' <C iV, thus accounting for the slope 1/T found 
on Fig. nil 

Turning to the case T < Tg/2, the asymptotic magne- 
tization can be computed from the following argument. 
Given the current state characterized by E and M, one 
can compute the probability that the magnetization takes 
the value M' after one transition, in the large N limit - 
see Eqs. and (|5^ in the Appendix. This distribution 

is gaussian and independent of E; its average value M 
can be identified with the most probable value: 



M-(l-§)M + 2i.(l-^l- (32) 



AT 



h 



The asymptotic magnetization can be computed self- 
consistent ly by looking for the fixed point M — M. 
Keeping only the first order in h, one can safely neglect 
the term M'^ /N'^, since one expects M ~ Nh. Solving 
the resulting equation and denoting the solution by M*, 
one finds M* = 2Nh/Tg. So, as here again C(t, i) = 1, 
the slope of the FDR is expected to be 2/Tg. 

So without entering into the details of the calculations 
developed in Appendix, this simple argument already al- 
lows to get an intuitive idea of the physical mechanisms 
responsible for the linearity of the FDR, and for the ef- 
fective temperature Tg/2. 



VI. DISCUSSION 

A. Effective temperature 

The effective temperature Tg/2 is very interesting for 
several reasons. On the one hand, it is surprising to see 
that this value is different from the usually expected value 
Teff = Tg which follows from mean-field spin-glass models 
|36j | and from Edwards- like arguments j3^ . To be more 
specific, one may expect from these theories a value 



T 



off 



dlnp 
dE 



(33) 



by associating lnp(_E) with the complexity (or configura- 
tional entropy) of the system. On the other hand, similar 
results have been found in a recent study of the Random 
Orthogonal Model (ROM) [s^, a spin-glass model with 
one-step replica symmetry breaking solution. Although 
the results on the ROM were not considered along the 
lines of the entropic-to-activated transition, we believe 
that it would be of great interest to search for such a 
transition in this model. To support this view, exponen- 
tial tails proportional to exp(— Q/A) in the distribution 
of heat exchanges Q have been found in the ROM ■ 
This suggests the existence of an underlying exponential 
distribution of energy levels with a tail: 



p{E) cx e^'^ 



E 



(34) 



An interpretation based on a scenario of spontaneous and 
stimulated relaxation in glassy systems , confirme d by 
numerical measurements of the FDR in the ROM |38l| . 
yields an effective temperature Toff related to A through 
Toff = A/2. In other words, the relation between the 
effective temperature and the slope of the exponential 
distribution is the same as in the NPP with 1 <C AT ^ 
N . This suggests that the mechanisms at play in the 
NPP should be rather generic, on condition that energy 
levels are exponentially distributed. However, this latter 
condition should not be so restrictive, since exponential 
distributions are expected for low lying energy states on 
the basis of extreme statistics [l3 • 



B. Influence of the energy distribution 

As already mentioned, two kinds of transitions between 
entropic and activated dynamics have to be distinguished 
in this model. On the one hand, a transition can be 
found as a function of temperature when crossing the 
value Tg/2. We have studied this transition in details 
within the NPP model. On the other hand, a crossover 
between both dynamics also appears, at fixed tempera- 
ture, when the energy of the system reaches the horizon 
level E^ . Here the regime changes as function of time 
rather than temperature. We did not study this phenom- 
ena in details within the NPP, as numerical simulations 
were difficult due to the large value of \E^\. 
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The former transition is controlled by the ratio ^ = 
T /Tg of the temperature T and the characteristic energy 
Tg of the exponential density of states p{E). On the 
contrary, the latter is due to the lack of direct paths be- 
tween states lying below the horizon level: the system 
has to reemerge first above the horizon, hence the onset 
of activated dynamics. So in the present model, differ- 
ent mechanisms are responsible for these two types of 
transitions. 

Yet, the first mechanism, which is related to the shape 
of the energy distribution, could also lead to a crossover 
as a function of time, for a given temperature. To this 
aim, one needs to consider a distribution p{E) which 
is not purely exponential. For instance, p{E) could be 
composed of a first exponential part p(E) oc e^/"^9 ' for 
E* < E < 0, and a second one p{E) oc e^^'^a^' for 

E < E*. Assuming T^'^'^ /2 < T < T^^^J_2, one then ob- 
serves an entropic dynamics as long as E{t) remains well 
above E* , and an activated one in the long time regime, 
when E{t) is below E* . 

Note that assuming on the contrary T^^^ < Tg'^\ one 
finds a temperature range Tg^^ 12 < T < Tg'^'' /2 such 
that the activated regime is found before the entropic 
one as the system ages (at constant temperature), which 
is rather counter-intuitive. Clearly, this example using 
an energy density with two different exponential parts is 
a bit artificial, but it is interesting for pedagogical pur- 
poses. A more natural example would be for instance a 
gaussian distribution: 

piE) = ^e-^'/'^o (35) 

The steady-state distribution at temperature T, P{E) cx 
p{E) e^^/^, is also a gaussian shifted towards low en- 
ergies, with average value Est = —Eq/T and the same 
variance E^. Assuming that we are in the low tempera- 
ture regime T <C -Eq, one has \Est \ ^ Eq. After a quench 
from high temperature to a low temperature T, the mean 
value E{t) is expected to drift slowly towards the steady 
state value Est during the aging regime, while the vari- 
ance remains essentially constant. When E{t) becomes 
close to Est -say E{t) ~ Est/2- one can linearize p{E) 
around E{t), to find locally an exponential distribution: 

P(£;,i)oce^/^(*) (36) 

with X{t) — EQ/\E{t)\. This exponential approximation 
is valid as long as \E — E{t)\ ^ Eq/T. In particular, 
if \E — E{t)\ « T, this approximation is fully justified. 
It is then natural to define a reduced parameter p{t) — 
T/X{t), which plays the same role as p in the BMM. 
Since \E{t)\ increases with time, so does p{t); equilibrium 
is reached for ^ = 1. When p{t) reaches the value i, 
one expects a crossover from an entropic dynamics to an 
activated one. 



This behavior should hold more generally for any dis- 
tribution p{E) of the form: 

p(£;) - e-"!^!" £;^-oo (37) 

with 7 > 1. On the contrary, distributions satisfying 
Eq. H37() with 7 < 1 should lead to the reverse behavior, 
that is an activated regime followed by an entropic one, 
at least in some temperature range. Note that such dis- 
tributions with 7 < 1 exhibit a glassy behavior for any 
temperature [l^ll5|. 



VII. CONCLUSION 

In this paper, we have established the dynamic phase 
diagram of the NPP model as a function of temperature 
and of the number K of spins allowed to flip simultane- 
ously. The first result is that some particular dynami- 
cal rules lead to the behaviors found in both usual trap 
models: the case K = 1 yields the activated trap model 
proposed by Bouchaud (BTM), whereas the model with 
K ~ N can be mapped onto the entropic version studied 
by Barrat and Mezard (BMM) . The former case has al- 
ready been studied in but was recalled here for the 
sake of completeness. 

In the intermediate range 1 <^ K ^ N , the dynamics 
of the energy remains essentially the same as for K — N, 
i.e. of BMM type, within the time window accessible 
with numerical simulations. For longer time scales, a 
horizon level Ef^ = —KTglnN is expected to yield a 
crossover to an activated regime, since states below Ej^ 
behave as isolated traps. Yet, the major difference with 
usual trap models is the presence of slow decorrelation of 
the observables: at each elementary transition between 
states, the magnetization decorrelates typically by a fac- 
tor {1 ~ K / N) , whereas full decorrelation in a single tran- 
sition is assumed for usual trap models. So the NPP can 
be mapped onto a modified version of the BMM which 
includes a slow decorrelation of the observable. 

This extra property has important consequences. On 
the one hand, in the entropic low temperature phase 
T < Tg/2, the correlation function decays only as [{t^j^ + 

t)/*^.]""^^^, with T] = {1- p)/{l-2p) and p = T/Tg, 
i.e. much more slowly than the hopping correlation func- 
tion. Indeed, in the thermodynamic limit TV — > 00, an 
ergodicity breaking occurs (the correlation function does 
not decay to 0) if K is such that K/N 0. Actually, 
for ^ 00, the correlation decays to at large time 
only if _ftr w aN, with a > 0. On the other hand, in 
the temperature range Tg/2 < T < Tg, the short time 
behavior of the correlation is precisely that of the usual 
BMM, with the onset of a singularity above Tg/2, but 
the long time (power law) tail becomes temperature de- 
pendent, with an exponent close to that of the BTM. 
This result suggests that thermal activation may be in 
that case the only relevant mechanism to decorrelate the 
observables, contrary to the activated phase of the BMM 
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in which both thermal activation and entropic slowing 
down control the system. 

In addition, we have studied the fluctuation-dissipation 
relation (FDR) in the NPP model. In the limit if < TV, 
the FDR becomes linear, and its slope depends on the 
temperature range considered. For T > Tg/2, the slope 
is equal to T~^, as in equilibrium and similarly to the 
case K = 1. On the contrary, for T < Tg/2, the slope 
is temperature independent and is equal to 2/Tg. Thus 
the effective temperature in this low temperature phase 

is Tg/2. 
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APPENDIX: Analytical developments on FDR 

In this appendix, we shall show, in the context of a 
modified BMM, how a well-defined effective tempera- 
ture Tg/2 emerges for slow decaying observables. To 
this end, we first study the observable 'magnetization' 
for a system microscopically composed of N spins 
and for which K spins at most can be flipped. Subse- 
quently, magnetization is independent of the energy and 
the energy evolution is the one given by the BMM en- 
dowed with Metropolis dynamics . We study the FDR 
for T < Tg/2, taking the limits in the following order: 
N —^ oo, /i — > and K oo. Then, in the same limits, 
we generalize the results to the general case of smooth 
observables that decorrelate by a factor (1 — K/N). 



FDR for the magnetization 

Zero temperature case 

We begin by studying the problem when temperature 
vanishes, choosing the magnetization as the (smooth) ob- 
servable. To derive the wished FDR, we need the rela- 
tion between the magnetization M before a jump and the 
magnetization M' after it. The system is assumed to be 
at an energy E, in the presence of an external field h. 

If K spins are chosen to be flipped with probability ^ , 
then K/2 spins are flipped on average. Assuming that 
K is large, fluctuations around this value K/2 can be 
neglected, and the effective dynamical rule is that K/2 
spins chosen at random are flipped (the new conflguration 
found in this way is then accepted or rejected according 
to the Metropolis rate). Given a magnetization M, the 



probability for a spin to be in the up state is given by 



1 M 

-2^2N 



(38) 



The probability to flip a number L < K/2 of up spins 
reads: 



(39) 



In the limit 1 <^ K <^ N, using the gaussian limit of 
Bernoulli processes, one finds for the new magnetization 
M' at each Monte Carlo step the following probability 
distribution: 



Pk{M'\M) 



1 / (Af'-(l-f)Af)2 
exp 



\/27raM 



2a 



M 



(40) 

with aM = SKpm{^—Pm)- Taking into account the zero- 
temperature acceptance rate, the distribution of magne- 
tization P{M'\E, M) after a transition is given by: 



PiM'\E,M) = pjsiMm. 



(41) 



r(£;, M, h) 

j dE' p[E') Q{E - hM - {E' - hM')) 



where is the Heavyside function, and T{E,M,h) is 
the normalization factor that can be interpreted as the 
escape rate from the state {E, M}: 



T{E,M,h) ^ J dM' pk{M'\M) X (42) 

J dE' p{E') Q{E - hM - [E' - hM')) 

One can then compute explicitly the distribution 
P{M'\E,M): 



P{M'\E,M) 



1 




(43) 



(Af'-(l-f)M-aM/3gM'^ 



2a 



M 



which, interestingly, appears to be independent of the 
energy E. This results in a decoupling between energy 
and magnetization. Note that apart from the normal- 
ization factor, the above calculation essentially amounts 
to multiply the a priori distribution pk{M'\M) by a fac- 
tor exp{f]gh{M' — M)). Given this distribution, we can 
recursively compute the magnetization of the system af- 
ter R jumps. Assuming that the magnetization is M^ 
at time when the field h is applied, and Mn after R 
jumps, one obtains the following relations: 



{MrM^)o 



a^{Ml), (44) 
a^(M^)o + 2p,Nh{l - a^) (45) 
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with a = 1 — K/N. The subscript indicates that the 
average is taken without any external field. Since there 
is no spontaneous magnetization, one has (Af^)o — 0. 
In this case, and in this case only, one finds using also 
{A'Oo = N: 



{MR)n = 2l3gh{{AOo - (MflAf„)o) 



(46) 



Strictly speaking, this relation is not the FDR since the 
parametric variable involved is the number R of jumps. 
In order to convert Eq. H46|) into a relation involving times 
tw and tw +t, we need to average it with the distribution 
of time intervals t given the number R of jumps. Let 
us consider this distribution in the presence of a small 
field h. At leading order in h, it is the sum of the zero 
field distribution and of terms proportional to Kh (we 
do not write the exact distribution here for the sake of 
simplicity) . We then see that the limits have to be taken 
in the following order: N oo, h 0, K ^ oo. Thus, 
the average of Eq. H46|l with that time distribution in the 
presence of h, combined with these very limits, leads to: 



distribution Pt{M'\E, M) now reads: 



Pt{M'\E,M) = 



PKiM'\M) X 



^E-h(M-M') 



(51) 



^jl^'^-~l3{E'~-hM'-(E-hM))^l3^E' \ 
I E~h(M-M') J 

where Tt^{E, M) is the escape rate at temperature T and 
with field /i, defined by normalizing P{M'\E,M). Per- 
forming the change of variables x — E' — h{M' — M) in 
the integrals, one gets: 

Pt{M'\E,M) = yj.^^'y ,^.HM'-M) , (52) 

( i-E MM'~M) \ 



d{M{t^ + t))h 



dh 



h=0 



TfT {{M{t 
^9 

(M {t^ + 



t)M{t^))n) 



Making the following identifications: 



C(i^,i^+t) ^ ^{M{ty,+t)M{t^))o 
we get the expected FDR, using also C{tw,tu]) = 1: 



(47) 

(48) 
(49) 



X{tw,t^+t) = —[l-C{t^,t^, + t)] (50) 

Note that the order of the limits does not restrict so 
much the validity domain of such a relation. Indeed, we 
have seen that K does not need to be very large to get a 
long BMM regime. And, for finite K, the same relation 
as Eq. can be exactly derivedj43. The reason we 
have chosen to present the large K solution lies in the 
simplicity of the relations between magnetization before 
and after a jump. In this case, it should be noticed also 
that this FDR is only valid for smooth observables with 
zero mean value at any time. Such prescriptions onto the 
observables are very similar to the ones needed to have a 
unique asymptotic FDR in the BTM [33 |. 



Finite temperature 

We shall see how to extend Eq. H5U|) to non-zero tem- 
peratures. In this case, Eq. (|40|l is still valid, whereas re- 
lation (|42|) is modified due to the Metropolis rates. The 



As a result, the dependence on h{M' — M) is essentially 
factorized, apart from the upper bound of the second in- 
tegral. In the regime T > Tg/2, this bound must play 
an important role, since the system visits quite often the 
high energy states. On the contrary, in the opposite tem- 
perature range T < Tg/2, these superficial states are no 
longer visited in the long time regime, so that one can 
neglect the infiuence of this upper bound and safely re- 
place it by since h < \M' - M\ ^ K. Within this 
approximation, the distribution Pt(M'\E, M) is again 
computed by multiplying pk{M'\M) by the exponential 
factor exp{f3ch{M' — M)), precisely as in the zero tem- 
perature case. 

Thus, for any temperature lower than Tg/2, the condi- 
tional distribution of magnetization after a transition is 
at large time the same as in the zero temperature case: 



Pt{M'\E,M) = P{M'\E,M) 



(53) 



Since moreover Pt{M'\E, M) does not depend on E, the 
derivation of the FDR is exactly the same as in the zero 
temperature case. This shows that an effective tempera- 
ture equal to Tg/2 is expected in all the entropic regime 

T < Tg/2. 



Generalization to other smooth observables 

More generally, one can consider observables A that 
decorrelate by a factor (1 — K/N) at each transition. 
In this case, we have the following relation between the 
value A of the observable before a jump and the value A' 
after: 



{^A), = ( 1 - f ) (A 



(54) 
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A natural stochastic evolution rule for such an observable 
is the following: 

^'=(i-^)^ + a (55) 

(^A is a gaussian white noise whose amplitude is imposed 
both by the conservation of variance of A, that is {A'^) — 
(A^), and by the correlation given in Eq. H54|l . Thus, at 
leading order in K/N, one can show that the variance of 
does not depend on the current value A. The variance 
is given by: 

(Ci) = ^(A^)o (56) 

Notice that we have implicitly chosen the case of observ- 
ables that are not correlated with the energy, since we 
assumed that (a does not depend on the energy. 



Then, we can see that Eq. H55() can be written as 
Eq. (gO)) with um = 2K{A^)/N. Thus, Eqs. gH) and 
(I45|l become: 

{ArA^)o = a^'iADo (57) 
{Ai^)h = a^(A„)o + 2/3,/i(A^)o(l-a^) (58) 

The end of the calculation is precisely the same as for 
the magnetization. 

To conclude, this analysis shows that a generic smooth 
observable is expected to satisfy a linear FDR with an ef- 
fective temperature Tg/2 when T < Tg/2. Moreover, we 
believe that the stochastic rule H55() is not an unavoidable 
ingredient for the emergence of the effective temperature. 
Indeed, as we already mentionned, il can be shown that 
magnetization verifies H5U|I even in the finite K case for 
which the gaussian law H40|l does not apply. 
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